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ABSTRACT 

66 ' 

^ . Observations of two H2CO (803 — 2o2 and 821 — 220) lines and continuum emis- 

sion at 1.8 mm towards Sgr B2(N) and Sgr B2(M) have been carried out with 
^ . the SMA. The mosaic maps of Sgr B2(N) and Sgr B2(M) in both continuum 

and lines show a complex distribution of dust and molecular gas in both clumps 
^ . and filaments surrounding the compact star formation cores. We have observed 

. a decelerating outflow originated from the Sgr B2(M) core, showing that both 

p |. the red-shifted and blue-shifted outflow components have a common terminal 

Q . velocity. This terminal velocity is 58±2 km s^^. It provides an excellent method 

^ in determination of the systematic velocity of the molecular cloud. The SMA 

'. observations have also shown that a large fraction of absorption against the two 

continuum cores is red-shifted with respect to the systematic velocities of Sgr 
, B2(N) and Sgr B2(M), respectively, suggesting that the majority of the dense 

00 ! molecular gas is flowing into the two major cores where massive stars have been 

^ ! formed. We have solved the radiative transfer in a multi-level system with LVG 

^ ! approximation. The observed H2CO line intensities and their ratios can be ade- 



O 



quately fitted with this model for the most of the gas components. However, the 



00 ! line intensities between the higher energy level transition H2CO (821 — 220) and 

^ ' the lower energy level transition H2CO (803 — 202) is reversed in the red-shifted 

outflow region of Sgr B2(M), suggesting the presence of inversion in population 
^ ! between the ground levels in the two K ladders (K_i= and 2). The possibility 

5^ I of weak maser processes for the H2CO emission in Sgr B2(M) is discussed. 

Subject headings: Galaxy:center — ISM:individual (Sgr B2) — ISM: kinematics 
and dynamics — ISM:molecules — radio lines: ISM — star:formation 
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1. Introduction 

The giant molecular cloud Sgr B2, located close to the Galactic center (~44 arcmin 
from Sgr A*), is a well-known massive star-forming region in our Galaxy. Sgr B2 consists 
of an extended envelope, a hot ring and a few compact cores (e.g. Goicoechea, Rodriguze- 
Fernandez & Cernicharo 2004). The radio continuum and recombination line observations 
of the compact HII regions suggest that Sgr B2(N) and Sgr B2 (M) are the two most active 
star forming cores in this region (Gaume & Claussen 1990; Gaume et al. 1995, Mehringer et 
al. 1993; de Free et al. 1995, 1996, 1998). Mascrs, outflows and possible rotation of the two 
dense cores have been revealed from observations of various molecular lines at centimeter 
and millimeter wavelengths (Reid et al. 1988; Gaume & Claussen 1990; Martin-Fintado et 
al. 1990; Mehringer, Goss & Falmer 1994; Lis et al. 1993; Kuan & Snyder 1996; Liu & 
Snyder 1999). In addition, previous observations have shown evidence for the two hot cores 
to be at different evolutionary stages and to have different molecular abundances (e.g. Vogel 
et al. 1987; Lis et al. 1993; Miao et al. 1995; Kuan, Mehringer & Snyder 1996; Liu & Snyder 
1999). 

H2CO pervades the interstellar medium and it has a simple chemical reaction path 
which has been proven to be a useful probe of physical conditions (e.g. Mangum & Wootten 
1993). The H2CO (lio — In) transition at 6 cm was observed in absorption against discrete 
continuum sources towards Sgr B2 complex with an angular resolution of ~ 10" x 20", 
showing nearly the same radial velocity pattern as that of the radio recombination lines 
(Martin-Fintado ct al. 1990; Mehringer, Falmer & Goss 1995). These authors suggested that 
the H2CO (lio — 111) transition probably arises from the surrounding gas with a relatively 
low mean H2 density of ~ 10"^ cm~^ (Martin-Fintado et al. 1990; Mehringer, Falmer & Goss 
1995). 

The millimeter H2CO lines are an excellent tracer of H2 density > 10^ cm~^ (e-g-, 
Mangum & Wootten 1993). In addition, H2CO is a planar asymmetric top molecule with 
very little asymmetry. The symmetry of the spin function of the molecule leads to two 
transition classes: ortho-H2CO levels if the spin wavefunction is symmetric and para-H2C0 
levels if antisymmetric. Since para-H2C0 is 1-3 times less abundant than ortho-H2CO, 
observations of para-H2C0 have less opacity effect (Kahane et al 1984; Mangum & Wootten 
1993). Hence, para-H2C0 appears to be a better probe to determine the physical conditions 
of the massive star formation regions. 

The millimeter/submillimeter transitions of H2CO gas require relatively high excitation 
temperature and high H2 density compared to those in the centimeter wavebands. If the 
brightness temperature of the continuum emission is higher than the excitation temperature, 
the absorption against the continuum cores can be observed in millimeter and sub- millimeter 
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wavebands with the high angular resolution of an interferometric array (such as the Submil- 



SMA, we have observed multiple H2CO lines towards Sgr B2 at 1.3 mm within a bandwidth 
of 2 GHz. Thus, with the same telescope system and calibration procedure, the uncertainties 
due to absolute flux density calibration among the different line transitions can be mitigated 
by measuring the line-intensity ratios, which are needed to determine physical conditions, 
such as kinetic temperature and H2 number density, of the gas. In addition, the SMA is 
not sensitive to extended larger scale emission (~ 50"). Thus, the SMA observations are 
sensitive to the clumps of high density gas rather than the extended diffuse components. 

In this paper, we present the results from the SMA observations of Sgr B2 at the 
H2CO lines and continuum at 1.3 mm. The paper is organized as follows: §2 discusses the 
observations and data reduction. In §3 we present the data analysis and results. In §4, we 
present the kinematics in Sgr B2(M) by a model incorporating a spherically symmetric inflow 
along with a decelerating outflow. In §5 we model the physical properties of the H2CO gas 
in Sgr B2 using the large velocity gradient (LVG) approach. In §6, we discuss the important 
results derived from our observations and analysis. We summarize the results in §7. We 
adopt a distance of 8 kpc to Sgr B2. 



Observations towards Sgr B2(N) and Sgr B2(M) were carried out at 218 (lower sideband) 
and 228 GHz (upper sideband) with seven antennas of the SMA in the compact-north array 
on August 1, 2005 for 8 hours. The projected baselines ranged from 5 to 50 kA. The 
weather was good during the observations with r~0.09 at 225 GHz. The typical system 
temperature was 130 K. Sgr B2(N) and Sgr B2(M) were observed in separate fields with 15 
min observing time on each source interleaving 5 min on the phase reference source Sgr A* 
(< 0.1 mas, ~3 Jy at 230 GHz, ~45 arcmin from Sgr B2). Our two target fields were centered 
at RA(J2000)=17M7°n9?882, DEC(J2000)=-28°22'18^'37 and RA(J2000)=17M7'°20?156, 
DEC(J2000)=-28°23'03r56, for Sgr B2(N) and Sgr B2(M), respectively. In addition, CaUisto 
(4.1 Jy) and the QSO 3C454.3 (32 Jy) were also observed for the flux-density and bandpass 
calibrations. The flux density was estimated from CaUisto with the assumption that its 
brightness temperature was 120 K and its angular size was 1''14. The three transitions 



"'^The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and 
the Academia Sinica. 




Taking advantage of the large bandwidth coverage of the 
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Observations & Data Reduction 
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H2CO (3o3 — 202), (322 — 221) and (821 — 220) were positioned in the lower 2 GHz sideband 
(LSB). H2CO (822 — 221) emission appears to be blended with a transition of CH3OH and will 
be not used in the analysis in this paper. The two unblended transitions, H2CO, (803 — 202) 
and (821 — 220), have upper level energies of 21.0 and 67.8 K and rest frequencies of 218.2222 
and 218.7601 GHz, respectively. The spectral resolution of 0.8125 MHz corresponds to a 
velocity resolution of 1.1 km s~^. 

The data reduction was carried out in Mznac?B We recomputed the Doppler veloc- 
ity for each of the target sources since the on-line Doppler tracking was only made on 
Sgr B2(N). System temperature corrections were applied. Antenna-based bandpass ripples 
were corrected by applying the linear interpolation of the bandpass solutions determined 
from Callisto and 8C454.8. There are spectral-window-based offsets in amplitude and phase 
on some baselines due to the correlator errors, which were also corrected. The residual errors 
due to the bandpass shape of the individual spectral windows were reduced to a level below 
1% of the continuum level. In the antenna-based gain corrections, we chose the visibilities 
of Sgr A* in the UV range between 20 kilo wavelengths and longer in order to eliminate the 
contamination from the extended dust and HII emission. The gains determined from the 
point source (emission from Sgr A*) were applied to the Sgr B2 data. 

The molecular lines were identified in the rest frame by use of the JPL catalog, as done 
by others (Sutton et al. 1985 and Nummelin et al. 1998). The H2CO transitions and a few 
other molecular lines have been unambiguously identified (see spectrum in Fig. 1). 

The continuum was subtracted with a linear fitting to the spectrum of line-free channels 
in each baseline. For Sgr B2(M), the line-free channels can be easily selected from the 
spectrum in the (u, v) domain (Fig. 1). However, for Sgr B2(N), the molecular spectral lines 
are crowded over each 2 GHz band and it is difficult to choose the line-free channels. We 
developed a procedure to select the line- free channels for the continuum subtraction. First, 
after reducing of the spectral resolution to 1 km s~^, we made each of the channel maps 
(2600 channels over the 2 GHz band) including both the line and continuum. From the 
channel maps, we selected the channels without extended emission (>4cr). The continuum 
level was determined from these apparently line-free channels. We selected 16 and 87 line- 
free spectral windows for the fields of Sgr B2(M) and Sgr B2(N), respectively. Then, using 
the task UVLIN in Miriad, we determined the continuum level by fitting to the line-free 
channels. The UVLIN gives two output (u, v) data sets, one for the continuum and the 
other for continuum-free spectral line. 



^The data reduction procedures for SMA data are outlined in the web site 
http: / /sma- www. cfa.harvard.edu/ miriadWWW, 
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Self-calibration was performed to the continuum data using the compact sources of the 
continuum cores for a few iterations in order to eliminate residual errors. The gain solutions 
determined from the continuum data were applied to the line data, from which we constructed 
the spectral data cube. The preliminary images of the continuum and lines were made using 
natural weighting. The clean algorithm was applied to remove the effects of the side-lobes. 
The synthesized beam sizes of the continuum and line images were approximately 5''4 x 3''2 
(PA=12.5°). The mosaic maps of the two fields, Sgr B2(N) and Sgr B2(M), were made using 
a simple linear mosaic algorithm. The primary beam attenuation was also corrected in the 
final line cubes and continuum image. The statistical la rms noise of the continuum was 
0.1 Jy beam~^. The la rms noise levels of the fine images were 0.18 and 0.16 Jy beam~^ 
per channel for the H2CO (3o3 — 202) and (821 — 220) lines, respectively. The higher noise 
in the H2CO (803 — 202) line images was due to its strong line intensity and the limit of the 
dynamic range in the clean process. 

3. Data Analysis &: Results 

3.1. The Continuum at 1.3 mm 

Fig. 2 shows the mosaic continuum map of Sgr B2(N) and Sgr B2(M) at 1.8 mm. The 
brightest components of Sgr B2(N) (Ip=29.2±2.1 Jy beam"!) and Sgr B2(M) (Ip=20.2±1.8 
Jy beam"^), are associated with the massive star forming cores Kl-3 and Fl-4 (Gaume & 
Claussen 1990), respectively. In addition to the emission from these cores, a few nearby con- 
tinuum clumps were detected, including the components K4 (Ip=1.40±0.12 Jy beam~^), NE 
(Ip=1.20±0.14 Jy beam-i), NW (Ip=1.00±0.11 Jy beam-^), ME (Ip=0.58±0.15 Jy beam-^) 
and MW (Ip=1.70±0.18 Jy beam-^). 

Gaussian fitting to the individual continuum components was carried out. The individ- 
ual emission clumps near the compact cores or the Kl-3 and Fl-4 clusters were modeled as 
simple Gaussian components. The two compact cores, Kl-3 and Fl-4, appeared to be too 
complicated to be fitted with single Gaussian components. A model consisting of a Gaussian 
and two unresolved compact components was used to fit to the data of the core (Kl-8). The 
peaks of the Gaussian component and one of the point components are consistent with K8 
and K2 positions, respectively. Another point component agrees (within ~ 1") with the 8.5 
mm continuum source that is located south of K3 (see Fig. 1(b) of Liu & Snyder 1999). Sgr 
B2(M) core (Fl-4) is fitted well with a Gaussian component and a point component (close 
to F3 within ~ O'.'S). The peak positions, deconvolved angular sizes, peak intensities and 
total fiux densities of the continuum components are summarized in Table 1. 
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The components K4, MW and Z10.24 have been detected at radio and millimeter wave- 
lengths (Lis et al. 1993; Kuan & Snyder 1994; Gaume et al. 1995; Liu & Snyder 1999). 
Z10.24 located in the middle between Sgr B2(N) and Sgr B2(M) (see Fig. 2) shows a unique 
filamentary structure at 1.3 mm continuum. The designation of Z10.24 follows that used 
by Gaume et al. (1995) and de Free et al. (1996) who detected the H66q; line towards it. 
Z10.24 was marginally detected at 1.3 mm by Lis et al. (1993) with higher angular resolu- 
tion (4'.'5 X 3'.'7) and poorer sensitivity. The SMA observations show an elongated structure 
in Z10.24 with a peak intensity of 2.02 ± 0.16 Jy beam^^ (12 a). The core of Z10.24 is 
unresolved in the sub-arcsec resolution images at both 1.3 cm and 3.5 mm (Gaume et al. 
1995; Liu Sz Snyder 1999). The 1.3 mm continuum image and the detected vibrational HC3N 
emission (de Vicente et al. 2000) suggest that the 1.3 mm continuum of Z10.24 is dominated 
by the dust emission and Z10.24 is likely to be a younger massive star formation region. 

There have been no detections of the components NW, NE and ME in the previous 
observations at longer wavelengths. These three components are possibly the dust emission 
from sub-cores at a relatively early stage of star formation. The detections need to be verified 
with further observations at shorter wavelengths and at higher angular resolutions. 

3.2. H2CO Lines 

The continuum-free channel maps in both H2CO (3o3 — 202) and H2CO (32i — 220) 
lines were constructed in the velocity range from 8 to 151 km s~^ at intervals of 1 km s~^. 
The channel maps of the H2CO transitions are complicated, containing several kinematical 
features in either emission or absorption. Those emission and absorption components are 
separated well in our higher spectral resolution maps but some of them (in the continuum 
core regions) are overlapped. In the moment analysis, the negative intensity value from the 
absorption and the positive value from the emission may cancel each other in the overlapping 
regions and thus the resultant moment maps might not reflect the true gas distribution. 
Hence, the emission and absorption need to be handled separately. 

3.2.1. Line Emission 

Fig. 3 is the integrated line emission images constructed from the channel maps for the 
two transitions, H2CO (32i — 220) and (3o3 — 202), respectively. The moment images were 
made with a 4 cr cutoff in each channel maps. The less signiflcant emission and the absorption 
are excluded in the moment analysis. Most of emission is distributed around the two cores of 
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Sgr B2(N) and Sgr B2(M). Clearly, the distribution of the H2CO emission is not spherically 
symmetric with respect to each of these cores. In the Sgr B2(N) region, in addition to the 
gas concentration at the core, gas clumps located north-east and south-west of the core are 
observed in both H2CO transitions. We note that in Sgr B2(N), the morphology of H2CO 
(821 — 220) emission, the higher transition gas, is similar to that of the emission from the 
lower transition gas H2CO (803 — 202). Hereafter we refer H2CO (821 — 220) as the higher 
transition and H2CO (893 — 202) as the lower transition. In the Sgr B2(M) region, a strong 
emission component elongated in northwest-southeast direction is observed in both H2CO 
transitions. An arch structure (Ml 30" long and 10" wide), 8" northwest of the Sgr B2(M) 
core appears in both H2CO transitions. The major difference in the emission distribution 
from the two transitions occurs in the outflow (Lis et al. 1998) region located south-east 
of the Sgr B2(M) core. A significant emission "tongue" (M5) (15" x 7") was detected from 
the higher H2CO transition gas while no significant detection was made of the lower H2CO 
transition gas. 



3.2.2. Systematic velocities 

Based on a line survey at 840 GHz from single dish observations, Sutton et al. (1991) 
obtained mean systematic velocities 65 and 61 km s"^ of Sgr B2(N) and Sgr B2(M). Their 
observations showed significant velocity variations among the different species. The differ- 
ences are mostly caused by the chemical differences of the molecules. The different species 
sample different physical environments. The high angular resolution observations of H66q! 
(de Pree et al. 1995; 1996) showed that the mean systematic velocities are 69.8 and 65.3 
km s~^ for the Sgr B2(M)-F and Sgr B2(N)-K clusters, respectively. The relatively higher 
mean velocity of 69.8 km in Sgr B2(M) is hkely caused by the high velocity motion of 
the ionized gas of the UCHII regions with respect to the centroid of the system. In section 
4, we will show that the systematic velocity determined from the terminal velocity of the 
outflow in Sgr B2(M) is 58 km s~^ (see also Appendix A). In the rest of the paper, we adopt 
systematic velocities of 58 km s~^ and 65 km s~^ for Sgr B2(M) and Sgr B2(N), respectively. 



3.2.3. Absorption towards the cores 

The strong continuum cores arc excellent probes of absorption by the cold gas that 
resides in front of them. Absorption is observed towards both Sgr B2(N) and Sgr B2(M) 
continuum cores (sec the top-left and bottom- left panels of Fig. 8). Multiple Gaussian line 
components were fitted to the spectra for each of the two transitions in both cores. The 
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parameters of these fits are listed in Table 2. The systematic velocities are marked with the 
vertical lines in Fig. 3. The majority of the absorbing gas in Sgr B2(N) and Sgr B2(M) 
is red-shifted with respect to the systematic velocities 58 and 65 km s~^, respectively. The 
red-shifted absorption gas provides evidence for the existence of gas accreting onto the two 
cores. 

If the absorbing gas covers the continuum source completely and the emission from the 
gas is insignificant, the line intensity is II ~ Ice~'^^. The optical depth (tl) can be derived 
from the formula 

rL = -Mr) = -Mi + ^), (1) 

ic ic 

where AIl = 1l — Ic is the observed line intensity and the Ic is the observed continuum 
intensity. The errors in Tl based on the fractional errors cil/Il of each line channel and 

(7ip/Ic if the line is optically thin ( — < 1 — ~ 0.6) are given by 

Ic 

a, ^ Vi^Tijhr + Micr. (2) 

In the optically thick case, the channels in the line center arc saturated and the linc-to- 
continuum ratio only gives a lower limit to the optical depth of the line. For example if 
n spectral channels are saturated by the absorption gas component at a velocity, we have 
II < 3(Til/ \/n in 3 a. Substituting this formula into the equation (1), the lower limit of the 
optical depth is 

V^Ic 

With the Rayleigh- Jeans approximation, the value of 1 Jy bearn"^ in our SMA observa- 
tions is equivalent to 1.5 K. The observed peak continuum intensities of the Sgr B2(N) and 
Sgr B2(M) cores are 29.2 and 20.2 Jy beam~^, which correspond to brightness temperature 
of 44 and 30 K in our observations, respectively. The solution of radiation transfer function 
in terms of the observed brightness temperature of the line (ATl^®) is given by 

^^obs ^ (j^Tex - foT'^^^l - e-'-), (4) 

where Tex is the excitation temperature of the molecular line and T^"^*^ = fc^c is the observed 
brightness temperature of the continuum emission with the true brightness temperature of 
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Tc; tl is the optical depth of the molecular cloud; for given solid angles of the molecular 

cloud (I^l)) continuum source (flc) and the telescope beam (fie), /l = + ^b) and 

fc = ^c/i^c + ^b) arc the beam filling factors of the line and continuum, respectively, if 
both the source and the telescope beam are in Gaussian shape; /o denotes the fraction of 
the continuum source covered by the molecular cloud. 

If the molecular cloud is in front of the larger continuum source (/o < "^), the observed 

J c 

brightness temperature of the line (Tl*^^) becomes: 

ATt = (ATe. - ^T°b^)(l - e-^). (5) 

JC 

If the molecular clouds cover the continuum cores with the same beam filling factor (/l = fc), 
the upper limits of the excitation temperatures H2CO absorbing gas ( ATl^^ < 0) are imposed 
by the observed continuum brightness temperatures of the continuum cores, i.e. 44//l and 
31//l K, for the Sgr B2(N) and Sgr B2(M) cores. For a given filling factor /l = 0.3, the 
excitation temperatures of H2CO would be less than 150 and 103 K for Sgr B2(N) and 
Sgr B2(M), respectively. 



3.2.4- Individual components 

Fig. 4 shows the spectra of the H2CO fines (averaged over one beam) made from the 
image cube for the rest of components. Each of spectral panels includes both the lower 
(3o3— 2o2, the green profile) and the higher (32i— 220, the red profile) transition lines. Gaussian 
fits to the H2CO spectra were carried out for both absorption and emission components. The 
emission and absorption probably come from the different regions along the line of sight. The 
angular resolution of our observations appears to be inadequate to distinguish the discrete 
components in the cores. However, our high spectral resolution is adequate to separate the 
emission from the absorption in the Gaussian fits. The parameters derived from the Gaussian 
fits to the two H2CO transitions are summarized in Table 2 including the central line velocity 
(Vlsr), the full width of half maximum (AV) and peak intensity (Ip). 

The Ml, M4 and M5 are located along the major axis of the bipolar outflow originated 
from the F-cluster (Lis et al. 1993). Ml is on the blue-shifted side of the outflow. Both the 
high- and low-transition spectra can be fltted with two Gaussian components at 52 km s"^ in 

emission and 64 km s~^ in absorption. The weak red-shifted absorption with respect to the 
systematic velocity suggests that a relatively cold gas component in front of the continuum 
source are moving towards it. The strong blue-shifted emission with respect to the systematic 
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velocity is the highly-excited gas emission in the outflow (in front of the continuum source) 
hkely mixed with the infaU gas emission (behind the continuum source). 

Sgr B2(M)-M4 is located close to the HII regions, i.e., the Fl-4 cluster. The higher 
transition spectrum can be fltted with two Gaussian components in emission at 66 and 70 
km s~^, both of which are red-shifted with respect to the systematic velocity 58 km s~^. The 
lower transition spectrum can be fitted with five Gaussian components at velocity 51, 65, 
70, 75 km s~^ in absorption and 100 km s~^ in emission. The morphology of this region is a 
complex. 

A possible model to interpret the spectral characteristics of M4 is considered here. 
If we assume a non-LTE condition for the gas and that both higher and lower transition 
gas comes from the same gas clump located in front of the continuum core, the excitation 
temperature of the higher transition gas is larger than the brightness temperature of the 
continuum, while the excitation temperature of the lower transition gas is less than the 
brightness temperature of the continuum. The red-shifted absorption suggests that the gas 
fiows towards the continuum core. The nature of absorption and emission at M4 suggests 
that the process of excitation of the infall molecular cloud is complicated. The radiative 
excitation by the strong FIR radiation field might play an important role in the region near 
the core since the collision alone can not produce the observed line ratio or the inverse 
population between the higher transition in K_i=2 and the lower transition in K_i=0 based 
on the LVG fitting (see Appendix B). 

The inversion in population between the two lowest K ladders is observed best in the 

red-shifted outflow component. At M5, the spectrum of the higher transition gas shows a 

very significant line emission (7-8 a), fitted to a Gaussian at 66 km s"^ with a line width 

of 16 km s~^, while less significant emission (1-2 a) of the lower transition gas is shown 

lf3 — 2 ) 

at the same position. The highly reversed line ratio jpr^ — t^H ~ 6 suggests that local 

thermodynamic equilibrium (LTE) is not valid at this location and a weak maser process is 
active in the outflow region. 

The gas components M2 and M3 and MW show signiflcant fine emission (> 10 cr) from 
the lower transition gas while the higher transition emission is relatively weak. The line ratio 

jhr^ — in those isolated component varies in the range between 1.5 to 5.2. 

H2CO maser at 6 cm (Mehringer, Goss & Palmer 1994) was detected in the Z10.24 re- 
gion. The SMA spectrum of the lower transition at Z10.24 shows that a signiflcant amount 
of gas is in absorption while the emission is present but red-shifted with respect to the 
absorption feature, a typical P-cygni proflle suggesting an outflow in this region. The spec- 
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trum of the higher H2CO transition can be fitted with two emission Gaussian components 
at 50 and 84 km with an absorbing Gaussian at 75 km s~^, which is consistent with 
the H66q; transition (de Prec, ct al. 1996). Excluding the possibihty of the expanding shell 
model, dc Prcc ct al. argued that an ionized outflow is likely centered at Z10.24. Our H2CO 
observations appear to favor their argument of a bipolar outflow from the UCHII region. 

Towards the NE continuum source, both the higher and lower H2CO transitions show 
that the majority of the gas in absorption is red-shifted with respect to the mean systematic 
velocity 65 km s~^. 

Towards NW, a broad (AV = 18 ± 1.3 km s~^) absorption (—1.1 Jy beam~^) from the 
lower transition line is detected at 71 km s~^. The spectrum of the higher transition shows 
no significant lines in either emission or absorption. 

3.3. Kinematics 

Figs. 5 and 6 show the images of intensity weighted velocity (or moment 1) of the H2CO 
emission gas from both transitions. The moment 1 maps were constructed with a cutoff of 
8 CT from each of the channel images in the velocity range of 8-151 km s~^. 

In the Sgr B2(M) region, the kinematical structure observed from the lower transition 
(see Fig. 5) consists of the highly red-shifted components 5" south-east of the compact 
core and a northeast-southwest arch structure. The component M4 appears to be a fast 
moving compact component (Vlsr ~ 100 km s^^) ejected from the core. The morphology 
of the northeast-southwest arch in Sgr B2(M) from the velocity field in the higher transition 
gas (see Fig. 6) appears to be consistent with that observed in the lower transitions. A 
velocity gradient is present south-east of M4 (shown in Fig. 6), which appears to indicate 
a decelerating outflow. In the case of Sgr B2(M), the higher H2CO transition appears to 
trace outflow well, which is consistent with the interpretation of the larger scale mass outflow 
based on the lower angular-resolution observations of NH3 and SO (Vogel, Genzel & Palmer 
1987). The higher angular- resolution observations of NH3 (Guame & Claussen 1990) showed 
that the red-shifted emission is located south of the F3 HII region, and the blue-shifted 
absorption is located north of the red-shifted emission, showing a velocity gradient in north- 
south direction. Based on their higher angular-resolution observations of OH maser and 
NH3, Guame & Claussen (1990) suggested that the north-south velocity gradient can be 
explained by a rotating disk or a torus of material with an extent of 2'' 5 surrounding the 
Sgr B2(M)-F cluster. Although the angular resolution in our observations is not adequate to 
verify the kinematical model proposed by Guame & Claussen (1990), our observed arch-liked 
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morphology of the H2CO gas shown in blue-shifted emission from northeast to southwest 
in the larger scale also indicates that the gas is undergoing a complicated infall process 
interacting with the outflow while the gas is spiraling onto the core rather than being in 
simple free fall. 

In the Sgr B2(N) region, the lower transition map (see Fig. 5) shows the red-shifted 
velocity located north of K2 and blue-shifted velocity located south of K2 with a north-south 
velocity gradient across K2. A similar velocity gradient has been observed in the HC3N line 
(Lis et al. 1993). These authors argued that the north-south velocity gradient in Sgr B2(N) 
traces rotation. The kinematics of the lower H2CO transition in emission gas of Sgr B2(N) is 
likely dominated by the gas rotating around the core. In addition to the north-south velocity 
gradient from the higher transition H2C0( 821 — 220) map as observed in the lower transition 
gas, an east- west velocity gradient is also observed in Sgr B2(N) (see Fig. 6). The outflow 
in Sgr B2(N) was observed in east- west direction (Lis et al. 1993). The higher angular- 
resolution observations of NH3 (Gaume & Claussen 1990) appeared to show a southeast- 
northwest velocity gradient across K2. The southeast-northwest velocity gradient in the 
larger scale from our observations appears to be consistent with the kinematical structure 
observed in NH3. Based on the higher angular-resolution observations of NH3, ruling out a 
simple outflow/rotation model, Gaume & Claussen (1990) suggested that several kinematic 
components in outflow, infall and rotation might be involved in the Sgr B2(N) core. The 
southeast-northwest velocity gradient observed from the H2CO emission gas in Sgr B2(N) 
appears to be caused by a combination of rotation, infall and outflow. The angular resolution 
of our observations is not adequate to discern the details of these motions. 



4. Outflow and Infall 

In order to better understand the on-going astrophysical processes in the star formation 
cores, we modeled the observed kinematics and the ratio of the line intensities. In comparison 
to Sgr B2(N), the kinematics observed from Sgr B2(M) appears to be relatively simple and 
characterized by infall and outflow. Fig. 7 shows the single fleld map of H2CO (32i — 220) 
centered at the F-cluster. The major axis of the red-shifted outflow (PA=158°) can be drawn 
by connecting the continuum core and the tips of the outflows (Ml and M5). 

Along the major axis of outflow (PA=158°), we have made a position- velocity diagram 
(PV) from the higher H2CO transition line cube (see Fig. 8). The solid contours show 
the emission and dashed contours indicate the absorption. This diagram shows two distinct 
decelerating outflow components clearly. The red-shifted component shows that the emission 
near the core (position at 0" in the vertical axis) starts with a high velocity of 106 km s~^ 
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and the velocity declines to a terminal velocity at 58 km s"^ as the gas goes to the outer 
region (20"away from the core). If wc shifted velocity to the terminal velocity 58 km or 
in the system rest frame, the decelerating red-shifted outflow is mainly located in the first 
quadrant with the absolute velocity decreasing from ~48 to km s^^. The third quadrant 
shows the decelerating blue-shifted outflow component with deceleration from ~50 to km 

Appendix A describes the numerical calculations that have been carried out to model the 
kinematical characteristics of a decelerating outflow combined with spherically infall. With a 
simple assumption of mass conservation and power law distributions in both outflow velocity 
and molecular density, the observed kinematics of the outflow components can be well fitted 
to the decelerating outflow model with an initial velocity of Vo=83 km s'^ at ro=l" (0.04 
pc) from the core center with the outflow opening angle of =30° and inclination of (p =45° 
for the red-shifted component (the thick red curve in Fig. 8). The blue-shifted component 
can be fitted with Vo=33 km s"^ at ro=l" (0.04 pc), 6 =30° and if =45° (the thick blue 
curve). The dashed curves range initial velocity Vq at ro=l" given the same opening angle 
and inclination angle. For the red-shifted outflow, Vq is from 250 km s~^ to 17 km s~^. The 
blue-shifted velocity ranges from 83 km s~^ to 4 km s"'^. 

Our observations and analysis have shown the presence of a decelerating outflow from 
the core in Sgr B2(M). Such a decelerating outflow can be caused by the cntrainmcnt of gas 
through the interaction between a fast wind flow from the central stars and an ambient cloud 
core, where the incorporation of additional mass into the outflow decelerates the wind (Lizano 
et al. 1988). As the wind moves farther from the cores, the outflow velocity terminates at 
the systematic velocity 58 km s~^ of the molecular cloud. We note that both the red- 
shifted and blue-shifted outflows terminate at a common velocity, which is close to the mean 
systematic velocity 61 km s"^ of Sgr B2(M) (Sutton et al. 1991). In fact, the common 
terminal velocity 58±2 km s~^ of blue-shifted and red-shifted outflow appears to provide 
an unambiguous method, independent of the chemical processes in molecular clouds, for 
accurately measuring the systematic velocity of Sgr B2(M). The uncertainty in the terminal 
velocity is mainly due to the velocity resolution in the PV diagram and the uncertainty in 
the model fltting process. 

Finally, the absorption located in the second quadrant is shown to be red-shifted with 
respect to the systematic velocity 58±2 km s^^. Thus, with the accurately determined 
systematic velocity we can be certain that this absorption feature arises from gas flowing 
towards the continuum source. 
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5. 



Excitation 



Based on the line ratio of the two transitions from the spectra, most of the gas in the Sgr 
B2 region does not satisfy the LTE condition. Under the assumption that the molecular cloud 
has spherical geometry, the non-LTE calculation requiring the collisional excitation rates for 
H2CO (Green 1991) has been carried out for Sgr B2 using the large velocity gradient (LVG) 
approximation (e.g. Evans, Davis & Plambeck 1979; Scoville & Solomon 1974; Lucy 1971). 
The radiative transfer calculations and model fitting are discussed in details in Appendix B. 
The fitting results for the line components for models with various sets of parameters are 
given in Table 3. 

We summarize the excitation properties for the various H2CO components in the fol- 
lowing: 

1. Towards the absorption components of the cores (Kl-3, Fl-4), the kinetic temperature 
depends on the beam filling factor of the continuum emission, ranging from a few tens K 
for a large continuum beam fc (0.3) to a few hundreds K for a small fc (0.05). High 
angular resolution observations are necessary to precisely determine the range of Tk. The 
typical values of are 10^~^ cm~^ for the core components. The derived column densities 
Npara-H2Co/'^V are in the range of 10^^~^^ cm~^ km~^ s. 

2. For the discrete H2CO components, no LVG solutions could be fitted to the observations 
with either model A or model B in which fc is large and the brightness temperature of 
the background radiation is too small to produce significant absorption lines as observed. 
For excitation temperature of a few tens K and higher, the intensity of absorption (see 
equation(5)) suggests that the beam filling factor must be small, i.e. fc ~0.1 or smaller, 
for a region with a continuum peak intensity of ~2 Jy beam~^ in our observations. The 
kinetic temperatures derived from the LVG fitting are all below 100 K for both emission 
and absorption components except for M3 emission component in which Tk = ISOlgQ^ K. 
In general the components outside the cores are cooler. The typical values of are 10^~^ 
cm~^ for the discrete components. The derived column density Npara-H2Co/'^V is in the 
range of 10^^~^^ cm~^ km~^ s. 

If3 — 2 ) 

3. For the components in the red-shifted outflow of Sgr B2(M), the line ratio _ 



is significantly smaller than unity, suggesting that the population of the ground levels in the 
two K ladders {K_i = and 2) is inverted. Hence, a weak maser process between the two 
different K ladders must occur in this region. Based on our radiative transfer code with the 
collisional excitation rates for kinetic temperature Tk< 300 K (Green 1991), we cannot rule 
out the possibility that the H2CO (821 — 220) line in Sgr B2(M) outflow is excited in the 
high temperature C-shocks. On the other hand, the strong FIR radiation field near the Sgr 
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B2(M) core might also play a role in the population inversion of the molecule. 

6. Discussion 

The assessment of molecular cloud mass from molecular lines can be affected by the 
excitation, opacity, abundance variations and gas dynamics of the molecular lines. The 
optically thin submillmeter dust continuum emission has been proven to be a good tracer of 
molecular cloud mass (Pierce-Price et al. 2000, Gordon 1995). If we take an average grain 
radius of 0.1 fim and grain density of 3 g cm~^ and a gas to dust ratio of 100 (Hildcbrand 
1983, Lis, Carlstrom & Keene 1991), the dusty cloud mass and column density are given by 
the formulae (Lis, Carlstrom & Keene 1991) 

where T is the mean dust temperature (K), Q[y) is grain emissivity at frequency z/, 5*,^ 
is the flux density corrected for free-free emission, O is the solid angle subtended by the 
source. Assuming Qiv) at 1.3 mm is 2x10"^ and the dust temperature is 150 K for Sgr B2 
( Carlstrom & Vogel 1989; Lis et al. 1993; Kuan, Mehringer & Snyder 1996), we derived the 
masses, H2 column densities and number densities. 

Because the continuum at 1.3 mm contains free-free emission, we estimate the phys- 
ical parameters using the flux densities of the continuum corrected for free-free emission. 
Assuming 3.6 cm continuum emission with FWHM beam of 3'.'8 x 2''0 (Mehringer et al. 
1993) of Sgr B2(N) (Kl-3) and Sgr B2(M) (Fl-4) come from optically thin free-free emission 
{Sy oc v~^-^\ we estimate the free-free contribution of 4.7 and 8.4 Jy (~ 9% and 24% of the 
total flux densities at 1.3 mm) at 1.3 mm towards the continuum cores Kl-3 and Fl-4 in 
our observations. Our estimates are consistent with the determinations of Lis et al. (1993) 
(6% and 33% of the total flux densities at 1.3 mm) and Martin-Pintado et al. (1990) (~ 
9% and 28% of the total flux densities at 1.3 mm ) for Kl-3 and Fl-4. The continuum flux 
densities at 1.3 mm corrected for free-free emission are 47.4 and 27.2 Jy for Kl-3 and Fl-4, 
respectively. 

From the flux densities at 1.3 cm (Gaumc ct al. 1995), the estimated free-free emission 
contributions at 1.3 mm are 0.02 Jy, 0.04 Jy and 1.26 Jy for K4, Z10.24 and MW, respectively. 
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The corresponding continuum flux densities corrected for free-free emission are 3.5, 9.1 and 
11.0 Jy, respectively. 

The derived H2 masses, column densities and number densities are summarized in Table 
4. The estimated H2 masses of the Sgr B2(N) core (Kl-3) and the Sgr B2(M) core (Fl-4) 
are larger than those given by Lis et al. (1993), while the H2 number densities are less than 
their results. This result is caused by the relatively larger size and higher flux densities of 
the continuum in our observations. 

The H2CO (3o3 — 202) spectra show absorption against both Sgr B2(N) and Sgr B2(M) 
compact cores and multiple absorbing peaks. The absorptions are dominated by red-shifted 
gas, suggesting that the lower transition H2CO (3o3 — 202) traces the cold gas in front of the 
continuum cores falling into the two compact cores. 

Previous observations showed multiple massive young stars in Sgr B2(N) and Sgr B2(M) 
(e.g., Gaume et al. 1995, de Pree, Goss & Gaume 1998). There are multiple absorbing peaks 
with different optical depths from our H2CO spectra in the Sgr B2(N) and Sgr B2(M) cores, 
which appears to indicate that the gas is falling into the massive stars or massive star forming 
cores embedded at different depths in the molecular clouds (Mehringcr, Palmer & Goss 1995). 
However, the angular resolution of our observations is inadequate for us to determine whether 
there are multiple regions present or whether the overall gravitation potential dominates the 
infalling gas. If infaUing gas is in simple free-fall, the infalling velocities can be estimated by 



where Tin is the infall radius, M is the sum of the gas and star masses included in the rj„ 
and G is the gravitation constant. 

The H2 masses derived from the continuum are 1.4x10^ and 7.9x10^ M© for the cores 
of Sgr B2(N) and Sgr B2(M), respectively. The VLA observations of radio continuum at 
1.3 cm (Gaume et al. 1995) showed that there are three UCHII regions (Kl, K2 and K3) 
in the core of Sgr B2(N) and four UCHII regions (Fl, F2, F3 and F4) in the core of Sgr 
B2(M). By use of the relationships between stellar spectral type and stellar mass (Vacca, 
Garmany & Shull 1996), a total stellar mass of 68 M© was inferred for the massive stars in 
the core of Sgr B2(N). The higher resolution observations (0''05) at 7 mm (de Pree, Goss & 
Gaume 1998) resolved out Fl, F2, F3 and F4 into twenty-one UCHII regions, and a stellar 
mass of 443 M© was inferred corresponding to the massive stars in the core of Sgr B2(M). 
The mass of the massive stars in the Sgr B2(M) core is six times larger than that in the Sgr 
B2(N) core. Taking the major axis sizes of 0.24 and 0.29 pc as the infall radii of the two 
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cores (Kl-3 and Fl-4), we inferred the infall velocities of 21 and 15 km for Sgr B2(N) 
and Sgr B2(M), respectively. Hence, based on our SMA observations, we have shown that 
high-density molecular gas is continuously feeding onto the active star formation cores in 
both Sgr B2(M) and Sgr B2(N). 

7. Summary 

The continuum emission and H2CO spectral lines at 1.3 mm were observed with the 
SMA. We detected the continuum emission from the compact cores of Sgr B2(N) and Sgr 
B2(M). Outside of the two compact cores, a few clumps, K4, NE, NW, ME, MW and the 
filament of Z10.24 were also observed. Except for the three newly observed objects, NE, NW 
and ME, all of the other compact cores, clumps and the filament have radio counterparts at 
centimeter wavelengths. 

The spectra in the two H2CO lines towards the continuum and outside of the continuum 
showed the gas components in either emission or absorption. From the integrated line in- 
tensity and intensity weighted velocity maps in line emission, we identified possible outfiow 
motion in Sgr B2(N) and Sgr B2(M) cores. The infalling gas was detected by the red-shifted 
absorbing gas against the continuum in Sgr B2(N) and Sgr B2(M) compact cores and the 
nearby dust clumps and filament Z10.24. The absorbing gas in Sgr B2(N) and Sgr B2(M) 
cores is red-shifted with respect to the systematic velocities, suggesting that the gas in front 
of the continuum cores is fiowing onto the continuum cores. 

In comparison to Sgr B2(N), the kinematics in Sgr B2(M) are relatively simple. We 
interpreted the observed outfiow and infall in Sgr B2(M) by a model incorporating a spheri- 
cally symmetric infiow with a decelerating outfiow. A decelerating outfiow from Sgr B2(M) 
was evident in our observations. We showed that both the red-shifted and blue-shifted out- 
flow components share a common terminal velocity. From the outflow terminal velocity, 
we determined the systematic velocity of 58±2 km s^"*^ for Sgr B2(M). With the systematic 
velocity well determined in a manner independent of chemistry in the molecular cloud, we 
are confldent that the majority of the absorbing gas at least in Sgr B2(M) is red-shifted and 
fiowing towards the active star forming core. 

Using the observed two H2CO lines, we have derived the excitation conditions of Sgr 
B2 region in LVG approximation. For the absorbing components in the Sgr B2(N) and 
Sgr B2(M) cores, the kinetic temperature, the H2 density and column density are in the 
range of several tens to hundreds K, 10^"^ cm""^, 10^^"^^ cm"^ km~^ s, respectively. The 
components outside of the cores are relatively cooler, and the H2 density and column density 
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are one order of magnitude less than those in the Sgr B2(N) and Sgr B2(M) cores. In the Sgr 

B2(M) outflow region, the intensity of the H2CO (821 — 220) is larger than that of the H2CO 
(3o3 — 202), suggesting the inversion of the distribution of the particle number between the 
two ground levels of K_i =0 and 2. The H2CO (821 — 220) might be excited by C-shocks in 
the outflow. The radiative excitation by the strong FIR radiation field in Sgr B2(M) may 
play an important role in the population inversion. 

We thank the SMA staff for making the observations possible. J.J. Wang acknowledges 
support from National Natural Science Foundation of China (No. 10328306). 

A. A Model for A Moleculeir Cloud With Inflow And Decelerating Outflow 

Components 

In this appendix, we discuss a model incorporating a spherically symmetric inflow (Shu 
1977) with a decelerating outflow (Cabrit & Bertout 1990, Raga et al. 1993). Wc calculated 
position-velocity (PV) diagrams based on this infall-outflow model. These PV diagrams 
can be compared with the observed PV diagrams from Sgr B2(M) to better understand the 
physical processes in the Sgr B2(M) molecular cloud. 

A.l. Infall 

For an ideal isothermal flow, under the assumption of spherically symmetric collapse, the 
volume density p(r) and velocity profiles V{r) as functions of radius r within the collapsing 
inner envelope can be described by the power-laws (Shu 1977), 

Pir) - Po(^)-S V{r) = -( ^)2, (Al) 

where R is outer boundary of the collapsing cloud in spherical geometry, po is the density at 
R, M[r) is the mass of the cloud interior to r, G is the gravitational constant. Because the 
actual core has a finite size, the infall is terminated at an inner radius of tq. The negative 
sign in the velocity indicates that the direction of the infiow is opposite that of the unit 
vector of radius in the spherical geometry. The geometry of a cloud with spherical infiow 
and decelerating outfiow is shown in the Cartesian coordinate system (x, y, z) in Fig. 9. 
The sky plane is the x-y plane. 
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The column density (N) at a sky position (x,y) can be calculated by integrating the 
volume density p{r) . The column density in the sky plane is 

N{x, y) = 2poR / (s^ + a'y-^ds, (A2) 

^0 



where s = z/R, a = y/{x/RY + {y/Ry. In order to get an analytical form for the column 
density, we express (s^ + a^)~4 in a Taylor expansion. Neglecting the contribution from 
higher order terms since the interesting region corresponds to a << 1 (or y^x^ + y'^ « R) 
in the practical case, the final solution of equation (A2) is 

p3/2 

N{x,y)^5po ^^,^y,^,^„ y^^T^O. (A3) 

At (x,y) = (0,0), the inflow terminates at the inner boundary Tq. Assuming the emission from 
the far-side gas is blocked by the central compact object, the effective column density along 
the hne-of-sight is 

7V(0,0)=2poi?[(^)-^-l]. (A4) 

Prom equation (Al), excluding the mass of the central compact object, the mass M(r) 
of the spherically symmetric cloud within r can be derived by 

M(r)= rAnr'pir)dr= T 4,rrVo(^r ^c^r = ^^^[(^)t - (^)i]. (A5) 

Jrn Jrn ^ 6 K K 



The mass interior to R is fti — 3"""^^ ^ "^^^^^f*^^® equation (A5) can be 

re-written as 



M(r)^M^(^)i. (A6) 

The radial velocity of infall gas is the projection of the velocity vector V on 2; axis. The 
radial velocity averaged along the line-of-sight is zero at a position away from the central 
position (x, y) = (0,0) due to the cancellation of the velocity in a spherically-symmetric inflow. 
At the position (x, y)=(0,0), ignoring the contribution of the emission from the region behind 
the central object, the velocity at (x, y)=(0,0) in front of the central object can be expressed 

by 
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_ jyir)dr _ JJ^^r^dr _ 4x/2 / gM^l - (j; 

Jro Jro ^ 



We note that the sign convention in the spherically symmetric description of infalling 
gas in front of the central object is opposite to the convention of the radial velocity with 
respect to the central object in the local standard rest (LSR) frame. In addition, in the 
derivation above, the systematic velocity is assumed to be zero. 



A. 2. Decelerating Outflow 

Molecular outflows are observed in star formation regions. Various models have been 
also proposed to interpret the observed outflow morphology (e.g., Raga et al. 1993; Cabrit 
& Berout, 1990). In this appendix, a PV diagram for a decelerating flow is calculated. The 
outflow velocity fleld and density are assumed to have power-law distributions along the 
major axis of the flow (cf. Cabrit & Bertout 1990) 

V{r) = VoC-^r, p{r) = poi^r, (A8) 
r r 

where po and Vq are the volume density and velocity at the inner radius r^j„, respectively; 
r = y/x^ + z^, the distance from the center, because the flow is in the x-z panel, i.e., y ~ 0. 
The geometry of the blue-shifted outflow is illustrated in Fig. 9 by the cone in ihex — z plane 
with an opening angle of 20. The inclination angle of (p is the angle between the z-axis and 
the axis of the blue-shifted outflow. The angle between the north and x-axis is the negative 
position angle of the blue-shifted outflow if we take the convention of position angle in AIPS. 
The radius of an outflow cross-section is equal to r ■tan(O). Given conservation of mass flux 
across a cross-section and < < 90°, the relationship of the power-law indices 5 and a 
can be derived 



6 = 2- a. (A9) 

The size of a cut across the major-axis of the outflow along the line-of-sight (2;-axis) can be 
expressed approximately by, 

Lz ^ tan(e)x/sin2((/7), (AlO) 
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if we take a — 1 and 5—1 (Cabrit & Bertout 1990), then the column density along the 
line-of-sight can be calculated from 



N{x) = J pdz = J po{^)dz, (All) 

the result of the integration is 

,w N N N , ,tan6 + (tan^O + sin^(;3)5 

sm (p 

The column density depends on the radius of the inner cross section (rmm)i volume density Pq 
at Tmin, the opening angle (0) and inclination angle (ip). In practice, the (/? and are in the 
ranges of < ip < 180° and < < 90°, respectively. A blue-shifted outflow corresponds to 
< 9? < 90° and 90° < </? < 180° to a red-shifted outflow. Fig. 10 plots g{Q, ip) as function 
of and (p. g{Q, (p) is a weak function of and (p. Given = 30° and (p — 45°, g is close 
to unity. 

The mean radial velocity of the outflow can be expressed as 

^ ^ J v{r)cos'pdz ^ ^r„,„J(0, ^) ^ g(0, v9)sinVcos(^ 

J dz x ' ' tan0 

The radial velocity of the outflow is inversely proportional to x along the major axis of the 
outflow and is a weak function (/(0, <^)) of and (p. As noted in the calculation of the infall 
velocity, the sign of the velocity used in the equation (A13) is opposite to the convention of 
the radial velocity with respect to the central source. The velocity of the central source with 
respect to the LSR, or systematic velocity, is not included here. 

Fig. 11 shows the plot of /(0, (p) as a, function of and cp. In the range of 90° < (p < 
180°, /(0, (p) has a negative value corresponding to a red-shifted outflow, while 0° < (p < 90° 
is for a blue-shifted outflow. Given = 30° and (p — 45°, / is about 0.6. In addition, for 
given (p, rmin, and Vq, the observed radial velocity decreases as increases. 

Assuming Vq = 85 km s"^ at Vq = 1" (0.04 pc), — 30° and (p — 45°, we can calculate a 
PV diagram of an outflow (the thick solid contour in Fig. 12). Fig. 12 shows a PV diagram 
for the parameter VoTminfi^, in a range between 0.5 to 5.1 km s~^ pc assuming that the 
blue and red-shifted outflows are symmetric around the central source, whose velocity is 58 
km s~^ with respect to the LSR. In this diagram, we also includes an absorption component 
of a spherical inflow with a mean infall velocity of 8 km s~^ with respect to the central 
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source. The PV-diagram (Fig. 12) shows the typical configuration of a spherically inflow 
and decelerating outflows in an active star formation region. 



B. Radiative Transfer with LVG approximation 



In order to model the excitation conditions and the physical properties of the molecular 
cloud components observed in Sgr B2, we solve for the radiation transfer in a multi-level 
system with the large velocity gradient (LVG) approximation. In this model, a molecular 
cloud is assumed to be spherically uniform at kinetic temperature (Tk), column density per 
unit velocity interval (Npara-H2Co/^V) and H2 density (nH2) with a volume fiUing factor /l y 
of the line emitting gas. The escape probability, /3 = (1 — e~'^)/r, is used to account for 
photon trapping. The volume filling factor /l,v is incorporated into the opacity calculation 
via the equation 



"° = W dVd^ ^"^^ - ^^'^ 

where tq is the optical depth at the line center, A21 is the Einstein coefficient for spontaneous 
emission, Xp;ira-H2CO is the abundance of para-H2C0, dv/dz is the velocity gradient, Xi and 
g[ are the fractional population and statistical weight of level i, respectively. The column 
density per unit velocity interval at the line center in the LVG model is given by 



IV /Av ^para-HaCPnHz/L.V .^^^x 
JMpara-HaCO/^V = dv/dz ' 

Given Uhj, T^ and Npara-H2Co/^V, the local LVG model gives the radiative temperature 
or the intensity of the line requiring as input the coUisional excitation rates. We adopted 
the coUision rates calculated by Green (1991). Green used the interaction potential between 
H2CO and He. Since He has twice mass of H2, coUisional excitation due to H2 should be 
2.2 times more effective than coUisional excitation by He (Green 1991; Mangum & Wootten 
1993). The coUisional excitation rates from Green were multiplied by this factor in our 
calculations. The accuracy in the Green's calculation of the total coUisional excitation rates 
is ~20%. 

Incorporating the Green's coUisional excitation rates and the volume filling factor dis- 
cussed above, we used the radiative transfer code in Miriad to solve for the radiative in- 
tensities from the para-H2C0 for the two transitions 803 — 2o2 and 821 — 22o- In the line 
emission case, the line intensity ratio of the two transitions is a good indicator of the ki- 
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netic temperature in molecular clouds {e.g. Mangum & Wootten 1993). This ratio becomes 
strongly dependent upon the volume density (uHa) when the kinetic temperature is greater 
than the upper state energy of the highest excitation transition (Eu). Given a column density 
(Npara-H2Co/^V), both the intensity of the lower transition I(3o3 — 2o2) and the intensity 
ratio (I(3o3 — 2o2)/I(32i — 220)) as a function of and can be calculated. We calculate 
five models with different sets of volume filling factor of the para-H2C0 gas (/l,v), beam 
filling factor of the line emission or absorption (/l) and beam filling factor of the continuum 
emission (/c)- In model A, B, C and D, the hue fiUing factors /l,v — 0.25 and ft, — 0.4 are 
assumed while /l,v = 1 and /l = 1 are assumed in model E. The continuum filling factors of 
fc — 0.3, 0.2, 0.1, 0.05 and 0.05 are assumed in model A, B, C, D and E respectively. Fig. 13 
shows the LVG results calculated with model C for the component Kl-3 at 72 km s~^. The 
absorption line intensities of the two transitions 3o3 — 2o2 and 821 — 220 are plotted as function 
of the kinetic temperature and H2 density. The line intensity ratio l(3o3 — 292) /I(32i — 220) is 
also shown (bottom panel). In the high density and cold region (the right-bottom corner), 
both transitions becomes optically thick and the excitation temperature is much smaller 
than the brightness temperature of the background continuum radiation and the line ratio 
becomes unity. When the gas becomes hot and less dense (left-top corner), the hue ratio 
becomes large. Fitting the observed results of the line intensity and the intensity ratio to 
the LVG model, we can determine T^ and uhj- We find that the absorption case is different 
from the case of emission, and the absorption line intensity of the lower energy transition, 
I(3o3 — 202), is a good indicator of the kinetic temperature (T^) of a cloud component while 
the line intensity ratio I(3o3 — 2o2)/I(32i — 220) places a strong constraint on the H2 volume 
density ( uhj) (see the top panel of Fig. 14). 

Fig. 14 shows examples of the LVG fitting with model C (/l,v = 0.25, /l = 0.4, and 
/l,v = 0.1) to the observed results from three typical regions. Top panel in Fig. 14 is 
the result showing the LVG model curves fitting to the observed I(3o3 — 202) and I(3o3 — 
2o2)/I(32i — 220) for the velocity component at 72 km s^^ in the region Kl-3 of Sgr B2(N). 
The dashed lines are absorption in the units of Jy beam"^. The thick dashed line and two 
thin dashed lines in light blue are the observed values of I(3o3 — 202) and 1 a uncertainty. 
The solid fines are the intensity ratio I(3o3 — 2o2)/I(32i — 220) ■ The thick fine and two 
tfiin lines in red correspond to the observed ratio and 1 a uncertainty. The hatched region 
marks the solution ranges of 2.9x10^-7.9x10^ cm~^ in nH2 and Tk = 350 — 360 K for given 
log(Npaia-H2Co/'^V) = 16.4 cm^^ km~^ s. We searched for LVG solutions in a wide range of 
11.0 - 19.0 cm-2 km-i s in log(Npa.a-H2Co/AV). 

The bottom-left panel in Fig. 14 shows the velocity component at 68 km s~^ in the Fl-4 
region of Sgr B2(M). Given log(Npara-H2Co/AV) = 17.0 cm^^ km~^ s, two possible solutions 
were found for this region, one with higher temperature (Tk — 150 K) and lower density 
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= 4.6 X 10^ cm ^ and other with lower temperature (Tk = 110 K) and higher density 
nH2 = 1.2 X 10^ cm-^ 

The bottom-right panel in Fig. 14 shows a LVG model to the observed results from the 
isolated emission region (M3) at 53 km s~^. Given log(Npara-H2Co/'^V) = 13.7 cm~^ km~^ 
s, Tk = 90 - 240 K and = 0.79 - 2.3 x 10^ cm~^ were found. The large uncertainties in 
Tk and are due to relatively weaker line emission. 
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Table 1. The SMA Measurements of Continuum Emission 



Source 


RA(J2000) 


ARA 


DEC(J2000) 


ADEC 


Angular size 


Peak Intensity 


Flux Density 




h m s 


(") 


/ // 


{") 


and (PA) 


(Jy beam~^) 


(Jy) 


Kl-3" 


17 47 19.89 


±0.1 


-28 22 17.6 


±0.3 


6;'8 X 3:'5(-25°) 


29.2±2 


52.1±4 


K4 


17 47 20.02 


±0.1 


-28 22 04.7 


±0.1 


5^2 X ^'.'8{-49°) 


1.4±0.1 


3.5±0.3 


NE 


17 47 20.56 


±0.5 


-28 22 14.4 


±0.3 


5"9 X i;'4(-71°) 


1.2±0.1 


2.5±0.3 


NW 


17 47 19.27 


±0.3 


-28 22 11.8 


±0.3 


5"1 X 4!' 1(0°) 


1.0±0.1 


2.3±0.2 


Z10.24 


17 47 20.04 


±0.3 


-28 22 41.2 


±0.4 


13f'3 X 4;'4(66°) 


2.0±0.2 


9.1±0.8 


Fl-4'' 


17 47 20.17 


±0.2 


-28 23 04.9 


±0.3 


7;'8 X 4;'4(-58°) 


20.2±1.3 


35.6±3 


ME 


17 47 21.65 


±0.6 


-28 22 57.0 


±0.7 


4!'l X 2;'6(73°) 


0.58±0.15 


1.0±0.3 


MW 


17 47 19.56 


±0.8 


-28 23 04.3 


±0.6 


12^8 X 8;'3(-13°) 


1.7±0.2 


12.3±1.3 



"Sgr B2(N) compact core (Kl-3) was fitted with a Gaussian source (6^8 x 3!'5, PA=-24.8°) centered on K3 
(RA(J2000)=17h47'°19?90, DEC(J2000)=-28°22'16:'9) and two point sources at K2 (RA(J2000)=17''47™19:'88, 
DEC(J2000)= -28°22'18;'6) and at (RA(J2000)=17''47™19!'95, DEC(J2000)= -28°22'12y3). The flux densities of 
the Gaussian source and the two point sources are 29.4±3.7, 19.1±1.5 and 3.6±1.0 Jy, respectively. The total flux 

density of the compact core (Kl-3) is 52.1±4 Jy. 

''Sgr B2(M) compact core (Fl-4) was fitted with a Gaussian source (7^8 x 4^4, PA=-58.1°) located 
within F cluster at (RA(J2000)=17'M7™20f 17, DEC(J2000)=-28°23'05;'l) and a point source centered on F3 
(RA(J2000)=17''47™20?14, DEC(J2000)=-28°23'04!'6). The flux densities of the Gaussian and the point sources 
are 21.3±2.6 and 14.3±1.0 Jy, respectively. The total flux density of the compa<;t core (Fl-4) is 35.6±3 Jy. 
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Table 2. The SMA Measurements of H2CO Lines 



H2CO (303 - 202) H2CO (321 - 220) 

Source Ipi Vlsri AVi n Ip2 Vlsr2 AV2 T2 

Jy beam~^ km km Jy beam~^ km km 



Kl-3 


5.9±0.1 


49±0.5 




-8.3±0.3 


56±0.5 




-19.0±1.6 


63±0.5 




-23.5±0.7 


72±0.5 




-20. Oil. 1 


84±0.5 




18.5±0.1 


96±0.5 




7.5±0.1 


106±0.5 




15.4±0.2 


121±0.5 


K4 


< 0.54" 




NE 


-1.2±0.1 


72±0.5 


NW 


-l.litO.l 


71±0.5 


Z10.24 


-l.OitO.l 


70±1.1 




1.2±0.1 


90±0.5 


Fl-4 


0.9±0.7 


45±2.0 




1.6±0.2 


54±0.5 




-6.1±0.4 


58±0.5 




-7.6±0.7 


62±0.5 




-16±0.1 


68±0.5 




-10±0.3 


76±0.5 




1.6±0.1 


96±0.5 




l.litO.l 


112±0.5 


ME 


< 0.76° 




MW 


6.6±0.1 


53±0.5 




-0.7±0.1 


70±0.9 


Ml 


6.7±0.4 


52±0.5 




-1.4±0.1 


64±2.3 


M2 


2.6±0.1 


55±0.5 


M3 


1.9±0.1 


53±0.5 


M4 


-2.1±0.1 


51±0.5 




-2.2±0.1 


65±0.7 




-0.8±0.3 


70±0.5 




-1.6±0.3 


75±0.5 




1.8±0.1 


100±0.5 


M5 


0.2±0.1 


66±3.1 



8±0.5 . . . 11.3±0.2 

3±0.5 0.3±0.1 4.7±0.1 

9±0.r5 1.1±0.2 -8.7±0.1 

12±1.2 1.6±0.1 -16.6±0.1 

11±0.5 l.litO.l -4.0±0.1 

7±0.5 ... 19.6±0.1 

8±0.5 ... 12.1±0.5 

15±0.5 . . . 7.8±0.2 

< 0.48" 

20±1.1 2..5±1.3 -0.5±0.1 

18±1.3 2.5±1.4 < 0.48" 

24±2.6 0.7±0.1 -0.3±0.1 

llitl.l ... 0.3±0.1 

18±7.8 ... 1.6±0.1 

2±0.5 ... 1.5±0.1 

3±0.5 0.4±0.1 -3.8±0.3 

4±0.5 0.5±0.1 -2.2±0.4 

10±0.5 1.6±0.1 -5.1±0.9 

3±0.5 0.7±0.1 -2.0±0.4 

17±0.9 . . . 1.7±0.1 

3±0.5 . . . l.OitO.l 

< 0.67" 

6±0.5 ... 4.3±0.1 

19±2.4 0.5±0.1 -O.liO.l 

11±0.5 ... 4.0±0.1 

21±3.4 0.9±0.1 -l.OitO.l 

6it0.5 . . . O.SitO.l 

7it0.5 ... 0.8it0.1 

7it0.5 O.SitO.l < 0.48" 

14itl.5 0.9it0.1 5.9it0.1 

3itl.l 0.2it0.1 4.2it0.2 

3it0.6 0.5it0.1 < 0.48" 

8it0.5 ... < 0.48" 

16it7.4 ... 1.2it0.1 



46it0.5 


18it0.5 




60it0.5 


5it0.5 




65it0.5 


4it0.5 


0.4±0.1 


73±0.5 


7±0.5 


0.8it0.1 


82±0.5 


4±0.5 


O.litO.l 


92it0.5 


7it0.5 




105it0.5 


lOitO.5 




116it0.5 


15it0.9 




71it0.9 


18it2.2 


0.5it0.2 


74it7.0 


15it9.5 


0.2it0.1 


87it6.5 


13it8.5 




39it0.5 


7it0.5 




53±0.5 


3it0.5 




58±0.5 


2it0.5 


0.2it0.1 


62it0.8 


4itl.4 


O.litO.l 


67itl.5 


5it3.4 


0.3it0.1 


74it0.6 


2itl.O 


O.litO.l 


91it0.5 


13it0.6 




107it0.5 


5it0.6 




53±0.5 


6±0.5 




70±3.8 


17±9.8 


O.litO.l 


52it0.5 


llitO.5 




64itl.2 


12it2.1 


o.aio.i 


53it0.5 


4itl.l 




55it0.5 


9it0.9 




66it0.5 


lOitO.5 




70it0.5 


3it0.5 





64it0.5 14it0.8 



" 3ct limit. 
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Table 3. LVG Model Results 



Source 


Vlsr Ipi/Ip2 Model" 


Tk 


1 / Npara-HoCO n 

log{ ) 




(km s-i) 


(K) 


(10^5 cm-3) (cm-2 km-l s) 



Kl-3 (a)*" 63 2.2±0.2 A ^^-^-tl 14.9-15.2 

B ISStg 17.8lJ°g'' 15.3-15.6 

C 200^^1° 10.0l^°5° 15.3-15.7 

D 300j:i5° 0.22t55;[5| 16.5-17.6 

E 255195 0-07lo:oi 17.1-17.3 

(a) 72 1.4±0.1 A 6ll^ 7ll^f 15.1-15.2 

B 105^35 8-31:1,1 15.8-16.2 

C 356tio 35t2?* 16.3-16.8 

D SOtll° 0.25tom 17.3-17.4 

E 1051^0 0-09io.M 17.3-17.4 

(a) 84 5±0.3 A 251^ 79l'23 13.7-14.1 

B 34l^ 316ligg* 13.8-13.9 

C 301" 316liJf 13.5-13.6 

D 200j:4° 0-025t55;og8 17.2-17.3 

E 200li°o O-OSio.oo 17.0-17.11 

NE (a) 72 2.4±0.5 C 25l:i 25.0ttl^ 14.6-14.8 

D 55t\o 6.3tg^'^ 14.8-15.8 

E 551^0 S.Ot^^g® 14.9-15.5 

Z10.24 (a) 70 3.3±1.1 C 26^1 lOtf 14.6-14.8 

D 50^20 '^■0-1%^ 14.8-15.6 

E 56li7 2.5t2'3'^ 14.7-15.1 

(e) 90 4.0±1.3 C 50tfl l-45io;6o 13.3-14.1 



D 



E 6OI15" 0.79tJ,;s9 13.0-13.9 

Fl-4 (a) 58 1.6±0.2 A 77I3 17.8155*5^ 15.4-15.7 

B I32I3 lOlf 15.7-16.0 

C 288t| 8.91^5^76 16.3-17.1 

D 575tJ° '^■OStill 17.1-17.4 

E 260t^° 0.05i!J;°i 17.05-17.15 

(a) 62 3.4±0.6 A 40^^^ 7.94+i^i.^| 14.6-14.9 

B 40^20 12.618^5* 14.4-14.7 

C ISOtgg 0.063t°;°43 16.88-17.04 

D 4:25tlf 0-28lo.?2 17.0-17.08 

E 280t?°(, 0.04lg;°J 16.86-17.02 

(a) 68 3.1±0.5 A 301?q '^^.Stfg^ 14.60-14.9 

B 41^11 31.6ti|* 14.5 14.7 

C llOlas 0.101:°;°^ 17.1-17.2 

D 465tf^5 O-l^lo'oT 17.0-17.18 

E 240tfQo O-OSto.oi 17.0-17.2 

(a) 76 5.0±1.0 A 26ljg 36.3^:^1^3 14.1-14.2 

B 25^12 63.112^3 13.8-13.92 



C 24IJ ISOI^g 13.25-13.76 

D 500tlo 0.018looo5^ 16.8-17.1 
E i50tro 0-02to:ooB 16.8-16.89 
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Table 3 — Continued 



Source 


Vlsr 


Ipl/Ip2 


Model" 


Tk 




log(— ' ) 




(km s-l) 






(K) 


(10^ cm-3) 


(cm~^ km~^ s) 


MW (e) 


53 


1.6±0.1 


C 
D 




4.50ll:« 


15.0-15.16 








E 


16+1 


7 94+7.86 


15.4-15.8 


(a) 


70 


3.5±1.7 


C 






14.83-15.77 








D 


52_22 




14.4-14.9 








E 






14.9-15.0 


Ml (e) 


52 


1.7±0.1 


C 
D 




o 47+1.54 
-1.07 


15.00-15.16 








E 


15+1 




15.3-16.8 


(a) 


64 


1.6±0.2 


C 


25+1 


■^•^'-1.07 


15.05-15.1 








D 


55+1 


20+1° 


15.2-16.1 








E 


90+1 


12.6+^1-?; 


15.2-16.0 


M2 (e) 


55 


5.2±1.0 


C 
D 


40il« 


io.o+_ii 


13.45-13.8 








E 


45l?g 


2.82l^;18 


13.00-13.85 


M3 (e) 


53 


2.4±0.3 


C 
D 




i-»O_0.66 


13.62-13.8 








E 




45+0 «l 


13.30-13.85 



"The filling factors are /c=0.3, /l=0.4, /l,v=0.25 for Model A, /c=0.2, /l=0.4, /l,v=0.25 
for Model B, /c=0.1, /l=0.4, /l,v=0.25 for Model C , /c=0.05, /l=0.4, /l,v=0.25 for Model D 
and /c=0.05, /l=1, /l,v=1 for Model E. For the absorption components, we used the following 
formula in the radiative transfer calculation; 

AT^''^ = (/LTex - |iT£b«)(l - e-^). 
JG 

For the emission components, we assumed /o = and used the following formula: 

ATf = = (/LTex - T,^b)(l - e-L), 



where T^mb is the temperature of the cosmic microwave background radiation. 
''The (a) and (e) indicate absorption and emission, respectively. 
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Table 4. The Properties of Sgr B2 Derived from Dust Emission 



Source 


Mh2 
(103 Mo) 


(1024 cm-2) 


(10^ cm-3) 


Kl-S" 


13.7±1 


21.2±2 


18.7±2 


K4 


l.OitO.l 


1.5±0.1 


1.3±0.1 


NE 


0.73±0.06 


3.2±0.3 


4.8±0.4 


NW 


0.66±0.06 


1.2±0.1 


l.libO.l 


Z10.24 


2.6 ±0.2 


1.7±0.1 


0.9±0.1 


Fl-4'' 


7.8±0.8 


9.0±1.3 


6.6±0.7 


ME 


0.28±0.03 


l.OitO.l 


2.4±0.3 


MW 


3.2±0.4 


l.litO.l 


0.46±0.06 



"The core of Sgr B2(N). 
''The core of Sgr B2(M). 
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217.5 218 218.5 219 

Rest Frequency (GHz) 

Fig. 1. — The spectrum from (u, v) domain on the basehne 3-4 for Sgr B2(M). The horizontal 
axis is the frequency in the rest frame and the vertical is the amplitude and phase. The 24 
spectral windows of the SMA's spectrometer are coded in different color. UL indicates the 
unidentified lines. We note that the features at 218.222 and 218.760 GHz are apparently 
dominated by H2CO (3o3 — 2o2 and 32i — 220) and the possible contaminations from O^^CS, 
C-C3H2 and SO^^O at the transitions near the H2CO frequencies 218.222 and 218.760 GHz 
appear to be insignificant in Sgr B2(M). These molecular lines also appear to be insignificant 
in Sgr B2(N). 
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17^47^22^ 21^ 20^ 19^ 18^ 

RA (J2000) 

Fig. 2. — The mosaic image of the continuum emission towards Sgr B2 at 1.3 mm obtained 
by combining the data from both sidebands (218 and 228 GHz). The synthesized beam is 
5'.'4 X 3'.'2, PA=12.5° (lower-right corner). The contours are -4, 4, 5.6, 8, 11, 16, 23, 32, 
45, 64, 91, 128, and 181 a. The rms (1 cr) noise level is 0.1 Jy bcam~^. The plus symbols 
indicate the positions of the continuum sources at 1.3 mm. The UCHII region Kl-3 and 
Fl-4 positions were taken from Gaume et al. 1995. 
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Fig. 3. — The mosaic images of integrated intensities with FWHM beam of 5'/4 x 3'/2, 
PA=12.5°. The pseudo-color is for integrated H2CO (32i — 220) emission intensity with 4 a 
cutoff in each channel (1 a is 0.16 Jy beam~^). The color scale on the right shows the range 
from to 80 Jy beam~^ km s"^. The contours (4, 8, 16, 32, 64, 128, 256 Jy beam~^ km s~^) 
are the integrated H2CO (3o3 — 202) emission intensity with 4 a cutoff in each channel (1 
a is 0.18 Jy beam~^). Top-left panel shows the observed spectra (solid curves) towards Sgr 
B2(N), green for H2CO (3o3 — 202) and red for H2CO (32i — 22o). The dashed curves are the 
Gaussian fitting to the multiple components. Bottom-left panel shows the spectra towards 
Sgr B2(M). The vertical lines mark the systematic velocities, 58 km s~^ for Sgr B2(M) and 
65 km s-i for Sgr B2(N). 
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Fig. 4. — Spectra from various components, green for H2CO (3o3 — 202) and red for H2CO 
(821 — 220) • Solid curves are the observed spectra and dashed are the multiple Gaussian 
fitting. The vertical lines mark the systematic velocities, 58 km s~^ for Ml, M2, M3, M4, 
M5 and MW in Sgr B2(M) and 65 km s^^ for NE, NW and K4 in Sgr B2(N). 



o 

17''47"'2lf5 21^0 20f5 20f0 19^5 19'0 18^5 18f0 17^5 

RA (J2000) 

Fig. 5. — The mosaic image of the intensity weighted velocity map in pseudo-color con- 
structed from the data after imposing a cutoff of 8 a, showing velocity field traced by H2CO 
(3o3 — 202) emission in pseudo-color. The scale on the right shows the velocity range from 
40 to 80 km s~^. The contours outline the integrated intensity of H2CO (821 — 220) emission 
at the level of 4 Jy beam"^ km s-\ The FWHM beam is 5"4 x 3'.'2 (PA=12.5°). 
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Fig. 6. — The mosaic image of the intensity weighted velocity map in pseudo-color con- 
structed from the data after imposing a cutoff of 8 a, showing velocity field traced by H2CO 
(821 — 220 ) emission in pseudo-color. The scale shows the velocity range from 40 to 80 km s^^. 
The contours outline the integrated intensity of H2CO (803 — 202) emission at the level of 4 
Jy beam-i km s'^ The FWHM beam is 5"4 x 8"2 (PA=12.5°). 
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Fig. 7. — The integrated line intensity of H2CO (821 — 220) in Sgr B2(M), showing the 
emission from the decelerating outflow along the direction in PA=158°. The color scale gives 
the intensity from to 90 Jy beam"^ km s'^. The FWHM beam is 5'.' A x 3''2 (PA=12.5°). 
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Fig. 8. — The velocity-position diagram cutting along the major axis of the outflows (see 
Fig. 7), which was constructed from the H2CO (32i — 220) spectral line cube. The solid 
contours show the emission from the decelerating outflows. The dashed contours are the 
absorption. The effective resolution of this map is labelled with the plus sign at bottom- 
right corner, i.e. 5 km s~^ in velocity and 5 arcsec in position. The vertical green line marks 
the systematic velocity. The solid curves are the best fitting to the decelerating outflow, 
with the decelerating outflow parameter Vormin/(0, = 2 km s~^ pc for the red-shifted 
component (red) and VQrminf{Q,(p) = 0.8 km s~^ pc for the blue-shifted. The dashed 
color lines show the range of Vor.mm/(0, f) between 0.4 to 6 km s~^ pc for the red-shifted 
component and 0.1 to 2 km s^^ pc for the blue-shifted component. 
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Fig. 9. — The geometry of the layout for a spherical infall (indicated by the small arrows 
towards the origin of the coordinates) and a decelerating outflow. Here the sky plane is 
the X — y plane. The outflow lies in the x — z plane with an opening angle of 2© and an 
inclination angle of </?. The z-axis points towards the observer. 




Fig. 10. — The geometric function g{Q,(p). 
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Fig. 11. — The geometric function f{Q,(p). 
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Fig. 12. — The PV diagram calculated from a model incorporating a spherical inflow with 
decelerating bipolar outflow for the parameter Vormm/(0, in a range between 0.5 and 5.1 
km pc (scaled in color). The position axis (vertical) is a cut along the major axis of the 
outflow. The velocity axis (horizontal) is the radial velocity of the outflow with respect to 
the LSR. The systematic velocity of the central object is assumed to be 58 km s~^ and the 
mean infall velocity of the absorption component is 8 km s~^. The contours correspond to 
Vormm/(6, V9)=0.8, 1, 1.4 and 2 km s~^ pc. The thick contour in the blue-shifted outflow 
corresponds to 2 km s~^ pc. The absorption (dashed contours) feature at V=66 km s~^ and 
X=0 arcsec shows the component of spherical inflow. The normalized absorption intensity 
is scaled in color. 
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Fig. 13. — An example shows the LVG analysis results derived with model C for component 
Kl-3 at 72 km s~^. Top panel: The line intensity of para-H2C0 transition I(3o3 — 2o2) as a 
function of and T^. Middle panel: The line intensity of para-H2C0 transition I(32i — 22o) 
as a function of and Tk. Bottom panel: The line ratio of I(3o3 — 2o2)/I(32i — 22o) as a 
function of and Tk. 



-44- 



K1-3 at 72 km s\ . 
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Fig. 14. — Examples of fitting the observed line intensity and intensity ratio to model C 
with the LVG approximation. Dashed lines are the line intensity of the Sqs — 2o2 transition. 
The solid lines are the hne ratio I(3o3 — 2o2)/I(32i — 22o)- The thick solid and dashed hnes 
correspond to the observed values of I(3o3 — 2o2) and I(3o3 — 2o2) /I(32i — 22o)- The hashed grey 
zones mark the solution regions in the UHa-Tk domain. Top: for the absorption component 
Kl-3 at 72 km s~^. Bottom-left: for the absorption component Fl-4 at 68 km s~^. Bottom- 
right: for the emission component M3 at 53 km s~^. 



